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Abstract 



We write down some relevant matrix elements for the scattering and decay 
processes of the pion by considering a quark-meson vertex function. The pion 
charge and transition form factors F n , i ? 7r7 , and i ? 7r7 * are extracted from these 
matrix elements using a relativistic quark model on the light-front. We found 
that, the form factors F n and F n ^ in the space- like region agree well with 
experiment. Furthermore, the branching ratios of all observed decay modes 
of the neutral pion, that are related to the form factors F^ and F n ^* in the 
time-like region, are all consistent with the data as well. Additionally, F n 
in the time-like region, which deals with the nonvalence contribution, is also 
discussed. 

PACS numbers: 12.39.Ki, 13.40.Gp 
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I. INTRODUCTION 



Form factors are important physical quantities in the understanding of the internal struc- 
ture of hadrons and relate to many other physical quantities. In this paper, we study two 
types of form factors for the pion: charge and transition ones. On the one hand, the former 
appears in the calculation of elastic electron-pion scattering in which one off-shell photon 
exchanges between electron and one of the quarks in pion. This form factor is related to 
the electromagnetic radius of the charged pion. The latter, on the other hand, comes from 
the reactions where the pion is coupled to two photons. This form factor also interrelates to 
the decay rates and branching ratios of all the observed decay modes of neutral pion. It is 
well known that these form factors must be treated with non-perturbative method. There 
are many different candidates for this purpose, such as lattice calculations vector me- 
son dominance (VMD) [@,[§, perturbative QCD (pQCD) with some non-perturbative input 
parameters 0-0], QCD sum rules P-pH|, and the light-front quark model (LFQM) |TT|-|T3|. 



LFQM is the only relativistic quark model in which a consistent and fully relativistic 
treatment of quark spins and the center-of-mass motion can be carried out. Thus it has 
been applied in the past to calculate various form factors 0-ff8|. This model has many 
advantages. For example, the light-front wave function is manifestly invariant under boost 
as it is expressed in terms of the momentum fraction variables (in "+" component) in an 
analogy to the parton distributions in the infinite momentum frame. Moreover, hadron spin 
can also be relativistically constructed by using the so-called Melosh rotation. The kinematic 
subgroup of the light-front formalism has the maximum number of interaction-free generators 
including the boost operator which describes the center-of-mass motion of the bound state 
(for a review of the light-front dynamics and light-front QCD, see ||19|| ). On the one hand, 
we will concentrate on the space-like region q 2 < (q being the momentum transfer) for 
charge form factors. In this region, the so-called Z graph contribution ||18[ vanishes and only 
the valence-quark contributes. As for the time-like region q 2 > for charge form factor, it 
concerns with the nonvalence contribution (Z graph). A reliable way of estimating this part 
is still lacking. On the other hand, the space-like and time-like regions, which correspond 
to the scattering and decay processes, respectively, are both considered for transition form 
factors. We make a consistent treatment with LFQM for the decay constants, the charge 
form factors, and the transition form factors which include the Z graph contribution within 
the time-like region. It must be emphasized that these derivations could be applied to all 
of the kinematically allowed region, no matter the momentum transfer is large or small. 



We compare these results with some experimental data on the charge form factors p0|-p3 
transition form factors |24 |, and branching ratios of some decay modes p5| for the pion. 



The paper is organized as follows. In Sec. II, the basic theoretical formalism is given and 
the decay constant, the charge and transition form factors are derived for the pion. In Sec. 
Ill, we fix the parameters appearing in the wave functions and calculate the form factors 
and branching ratios. Finally, a conclusion is given in Sec. IV. 
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II. FRAMEWORK 



In this section, we first write down some relevant matrix elements in the covariant form. 
After integrating out the minus component of the internal monentum, we introduce the 
light-front quark model to formulate the decay constant, charge form factors and transition 
form factors. 



A. decay constant 

The decay constant of a pion is defined by 

(0|A"|7r(p)) = V2iUp>>, 



(2.1) 



where = ?2T M T5?i is the axial vector current. Assuming a constant vertex function 
|jnyi6(| which is related to the qq bound state of the pion. Then the quark-meson diagram, 
depicted in Fig.l (a), yields 



(0 1 927^75?! k(p)> 



N c A n 



Tr 



75 



ivp\— i> + m) 



(Pi ~~ P) 2 — rn 2 + it Pi — rn 2 + ie 



7 M 75 , 



(2.2) 



where m = m u = rrid and N c is the number of colors. In terms of the LF coordinates 
(Pi,Pi,Pi±), we obtain 



|g27 M 7sgik(p)) 



X 



2(2tt) 4 
[Pi -p 



(4p^ 



dp i dpid 2 pi ± 
Pi(Pi-p) + 
m 2 + (pi — p)' 



/ _ m 2 + pf ± ie N 

Pi + I + 

v Pi Pi j 



+ 



ie 



n -l 



(pf-p+) (pt-p + ) 



(2.3) 



Comparing ( |2.1| ) with ( |2.3| ) and performing the integration over the LF "energy" p 1 in ( |2.3|) , 
we obtain 



dp\d 2 p\ 



rn 



2(2tt)3 pt(p + -pt) 



K[p~ 



m 2 + (pi — p)\ rn 2 + p\ 



P ~ Pi 



Pi 



(2.4) 



B. charge form factor 

The charge form factor of the pion in the space-like region, as being illustrated in Fig.l 
(b), is determined by the scattering of one virtual photon and one pion. This form factor 
can be defined by the matrix element 

(n'(p')\J^(p)) = eF:(q 2 )(p + pr, (2.5) 

where J M = qe q e^q is the vector current, e q is the charge of quark q in unit of e, q 2 = 
(p' — p) 2 < 0, and the superscript "s" represents the space-like region. From Fig.l (b), we 
obtain 
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i{p° 



(27 

h + m) 



(p' — pi) 2 — m 2 + it ^ (p — pi) 2 — m 2 + it 
Using the LF coordinates and taking the trace, we obtain 

eA^A^ f dpidpid 2 pi ± 



i(— p*\ + m) 

TT— 75 

Pi — m z + it 

iip 1 — i>\ + m) 



(2.6) 



(tt (p)|ge g e7 M g|7r(p)) 



2(2tt) 4 J p+(p- pi )+(p>- pi )+ 



xlp -p l 



m 2 + (p — pi) 5 



J" 



7 6 



m 2 + p 2 ± it 



Pi 



Pi 



P ~ Pi P ~ Pi 

2 i f^l ^ \2 



x p'-- ? r- - + , + (p 'T )x + 



n -1- 



where 



P = 4 



p /+ — 



Pi (Pi — P ■ p' — m 2 ) + P^(—Pi + Pi ■ p' + rn 2 ) + p'^(—p\ +p% ■ p + m 



(2.7) 



(2.8) 



Let us extract F n (Q 2 ) from the "+" component of the vector current J M . It does not loss 
the generality if q + is set to be zero for the momentum transfer in the space-like region. 
Performing the ^-integration, we obtain 



(e u + eg)I 



dpjd 2 pi L j 

2(2tt)3 \pt(p-Pi)+(pf-pi) 



K\p~ 



m 2 + (p — P\)\ m 2 + p\ 



P ~ Pi 



Pi 



x 



K>\p' 



m 2 + (p' — pi)\ m 2 + pl ± 



P l+ ~ Pi 



Pi 



(2.9) 



where 



1+ 



2p^ 



q+=0, p~ 



(2.10) 



By contrast, when we consider the process 7* — > tttt which is illustrated in Fig. 1(c), 
the momentum transfer q 2 > AM 2 is in the time-like region. In this region, the charge form 
factor is determined by another matrix element 



(n'{p')n{p)\qe^q\Q) = -(e u + e^eA^A^ 



d 4 pi 

(2tt)' 



Tr 



+ m) 



x 



i(p / — i>\ + m) 



Pi — m z + tt 
i{4— tfi + m 



75 



(p _ Pi) 2 — m 2 + ze (g — pi) 2 — m 2 + it 



(2.11) 



where q = (p + p'). It is well known that q 2 = q + q~ — q\ in the LF coordinate, if we firstly 
assume (^TJ the momentum transfer is purely longitudinal, i.e., q± = 0, q 2 = q + q~ > will 
be ensured. The same results will be obtained if one carries out the integral and then takes 
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q± = 0. However, the former method will reduce the calculating processes considerably. 
Thus this form factor is obtained as 



d 2 Pi ± [P 



m 2 + {q — p 1 ) 2 L m 2 +pf 



2(2vr)3 7 yi pt{q-Px) + {p-pi)- 



x 



K\p~ 



m 2 + (p — Pi)\ m 2 + p1 ± 



q+ -pi 
-l- 



p + — pi 



pi 



d 2 Pi_ 



dpi 



KoA 



2^7T 



X 



m 2 + (q-pi) 
q+ -pi 

,_ m 2 + (p-pij]_ m 2 +p 2 ± \ h 



pi 



m + pi ± \ 



2(2tt) 3 V ^ pi(q- Pl )+( Pl -p)+ 



Pi 



pi-p+ (q-pi)' 
where the superscript "t" represents the time-like region, 



K+ 



2g^ 



K+ 



Pi = — + 



and 



= 4 



Pi(m 2 — p\ + 2p\ ■ p — p ■ q) + p^(—p 2 i + pi ■ q + m 2 ) + q^(pf — p\ ■ p — m 2 ) 



L' K2= W 



(2.12) 



(2.13) 



(g-pi)~ 



(<j-pi)+ 



. (2.14) 



C. transition form factor 

There are two types of transition form factors, F nl and F nl *. The former comes from a 
neutral pion which coupled to two photons with one on-shell and the other off-shell, while 
the latter arises from the 7r°7*7* vertex, where 7* represent an off-shell photons. These 
two form factors are related to many interesting experimental data. In this section, we will 
calculate them in detail. 



1. transition form factor F n 



7TJ 



The form factor F ny is defined by the 7177* vertex 



-ie F nJ (q )e iivpa p v q p e° , 



(2.15) 

where q(e) is the momentum (polarization) of the on-shell photon. This process is illustrated 
in Fig.l (e). The amplitude is given by 

r M = -e (? e 5 /e 2 y / iVcA 7r J ■ 1 1 



x|Tr 
+ Tr 



75 



(2tt)< 
+ m) 



i(p' 1 - p 1 - 4+m) 



iip'i + m) 



(pi — p) 2 — rn 2 + it (pi — P — q) 2 — rn 2 + Pi — m2 + 



75 



i{A+ 4+ m ) 



(Pi — P) 2 — m 2 + ie (pi + q) 2 — m 2 + ie p\ — m 2 + ie 



iip'i + m ) 



(2.16) 
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Using the light-front coordinates and taking the trace in( 2.16|) , we obtain 

dpidpfd 2 pi ± 



T M = e (? eg/e 2 y / iVcA 7r y ■ 



2(2tt)< 



-Mme tlvpa tfq p E a 



pt(pi - p) + (pi + q) + \Pi -v 



m 2 + (pi — p) 



± 



if 



x 



+ 



X 



Pi 



m 2 + p1 ± if 



Pi 



Pi 



Pi +(f 



pi -p+ 
m 2 + (pi + q)\ 
Pi + Q + 



+ 



if 



pT + <r 



ptipi - p) + {pi - p - q) + \Pi - p - 



f _ m 2 + p\ | if \ ( 

Pi + — + — \\Pi - p -q 

V Pl PlJ\ 



m 2 + (p\ — p)" 
Pi ~~ P + 



+ 



if 



Pi -pi 
m 2 + (pi - p - q) 2 ± 



if 



Pi — p + — q + 



Pi — p — q 



(2.17) 



F W1 can be extracted easily by comparing ( |2.15| ) with fl2.17| ). When we consider q 2 < 0, the 
^-integration is performed in the condition q + = and the result is 



F^(q 2 ) = -A^N c e q e,, J 



dpfd 2 pi ± 
2(2vr) 3 ' 



Ajp~ 



m 2 + (pi — p) \ m 2 +p 2 ± \ 1 



x 



+ 



pt 2 (p-Pi) + [(p + q) - 
pt(p-pi) +2 \(p + q)~ - 



m 2 + (pi-p) 2 ± m 2 + (pi + q) 2 ± 
p + _ p + pf , 

m 2 + (pi — p — q) 2 ± m 2 + (pi + q) 



P ~ Pi 



Pi 



(2.18) 



If the decay process q 2 > is considered, the momentum transfer is in the time-like 
region. The pj~-integration is performed in the condition q± = in two regions: < pi < q + 
and q + < pi < p + , and in their photons-exchanged part: < pf < (p — q) + and (p — q) + < 
Pi < P + 



F'(q') = AJN c e q e^m 



d 2 pi ± 
2(2tt) 3 '1 Jo 



x 



2 i ^2 



m 2 + (pi — p)x m 2 +p 2 ± \ 1 



p^ -p{ 

2 , ^2 \ T -1 



Pl 



+ ( \ + ( N+/ - ™> +Pl± m + Pjx 

Piip-pi) yq-Px) \q + - — — 

V Pi q + - Pi 



+ / dpt 



K\p - 



m 2 + (pi — p) 2 ± m 2 +pf ± \ 1 



p^ -p{ 



Pi 



x 



+ 



X 



pt(p - pi) + (pi - q) + (p - q)' 



m + (p\ — p) j_ m + pf ± 



P + ~ Pi 



Pi ~ T 



v -<?' 



dpi 



K\p~ 



pf{p - pi) i (p - q - pi) (p - q) 



m 2 + (pi — p)\ m 2 + pf ± \ 1-1 
p+ -pl pi ) 

m 2 + (p — pi) 2 ± m 2 + (p — piY 



{p-piY 



(p-q-pi)- 
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+ / dpt 

'V q. 



Att \P~ 



m 2 + (j>i — p)\ m 2 + p\ 



P 



pf(p-pi) + (pi-p + q) + \q - 



Pi Pi J 
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(pi - p + q)~ 



(p-piY 



(2.19) 



transition form factor F n ^* 



As for the form factor -F^*, it is defined by the ^7*7* vertex |26 

r M „ = -ze 2 ^.^ 2 ,^ 2 )^^^, 



(2.20) 



where Q = ^(q' — q) , V = q' + q, and q' = p + q. This process is depicted in Fig.l (f). The 
amplitude can be written as 



-e q e q ie 2 \l N f A % 



x|Tr 
+ Tr 



75 



4+711) 



{pi — p) 2 — fn 2 + (pi — P — q) 2 — m 2 + ie lv p\ — m 2 + ie 



H01- 



m) i(pi + m) 



75 



p 1 + m) 



"7 



Kp'i+ (/ + m ) 



7m; 



z(^i + m) 



(pi — p) 2 — m 2 + ^ (pi + <?) 2 — m 2 + it Pi~ m2 + 



(2-21) 



The derivation of F^* is almost the same as that in ( |2.18||rT9"|) . But there is one thing 
needed to be clarified: (p + q)~ in (|2.18|J2.19| ) is the LF energy component of the on-shell 
photon, while it is the one corresponding to the off-shell photons in F%L. This difference 
will become more clear in subsection D. 



D. Light-Front Formalism 

Within the light-front formalism, a meson consisting of a quark qi and an antiquark g 2 
with total momentum p and spin S can be written as 



\M(p, S, S z )) = J {d 3 Pl }{d 3 p 2 } 2(2tt) 3 5 3 (p -fa- P2) 



E * S *"(Pi,P2,Ai,A 2 ) |gi(pi,Ai)g 2 (p 2 ,A 2 )>, (2.22) 



where p\ and p 2 are the on-shell light-front momenta, 

,2 j „,2 
P 

and 



P = (p + ,P±) , P± = (P 1 ,P 2 ) , P~ = (2.23) 
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1 Pl - 2(2vr)3 ' 

\q(pi, Ai)^, A 2 )> = (Pi)4 2 (P2) |0), (2.24) 
{&A'(P'),4(P)} = {<W),4(p)} = 2(2tt) 3 <5 3 (p'"P) * va . 
In terms of the light-front relative momentum variables (x, k±) defined by 

Pi — (1 — x)p + , P2 = xp + , 

Pi± = (1 - z)p± + k±, p 2 ±_ = xp±_ - kx, (2.25) 

the momentum-space wave-function ty SSz can be expressed as 

* 55 «(pi,p 2 , A 1; A 2 ) = (a, 0(x, fcj, (2.26) 

where <f>(x, k±) describes the momentum distribution of the constituents in the bound state, 
and R\A„ constructs a state of definite spin (S, S z ) out of light-front helicity (Ai, A 2 ) eigen- 
states. Explicitly, 

R S x %(x,k ± ) = ^(A 1 |^ / (l-x,fc ± ,m 1 )| Sl )(A 2 |^[ / (x,-A ;± ,m 2 )| S2 )(i Sl ;i S2 |5,S z ), (2.27) 

si,s 2 1 1 

where |sj) are the usual Pauli spinors, and TZm is the Melosh transformation operator: 

v t 1 \ m i + xM o + ia-k±xn 

Km(x, k±, nii) = (2.28) 

(rrn + xM ) 2 + k\ 



with n = (0, 0, 1), a unit vector in the ^-direction, and 



Ml = -± ± + 2 ± . 2.29 

(1 — x) x 

In practice it is more convenient to use the covariant form for R\,\ 2 [14[: 

R S x?x 2 (x, k ± ) = — u( Pl , Xi)Tv(p 2 , A 2 ), (2.30) 

V 2 Mq 

where M = ^Mq — (m x — m 2 ) 2 , T = 75 (pseudoscalar, S=0), and 
^ u(p, X)u(p, X) = p' + m, u(p, X)u(p, A') = 2 mt^y , 

A 

J^u(p, X)v(p, A) = jrf- m, v(p, X)v(p, A') = —2 m5x,\>. (2.31) 

A 

We normalize the meson state as 

(M(p', S', S' z )\M(p, S, S z )) = 2{2-nfp + 5\p' - p)5 s , s 5 s , Sz , (2.32) 

so that 



dx d 2 k_i 
2(2vr) 3 



(x,k ± )\' 



(2.33) 



After considering the normalization condition, the momentum distribution function 
(x, k±), for the pion, is related to the bound state vertex function A n by fl4| , |16 



A, 



(P 2 - Ml 



x(l — x) 
s/2 M 



-<j> v (x, k±). 



(2.34) 



Now, we can use the definitions of the light-front relative momentum variables (x, k±) ( |2.25 ) 
and the relation between the constant vertex function A n and the momentum distribution 
function <p(x, k±) (|2.34 ) to rewrite the equations of those decay constants and form factors 
in the subsections A, B, and C. The results are the following. For the decay constant, 



f* 



2\ N, 



dx d 2 k 



2(2tt) 3 

For charge form factor in the space-like region, 



<j> v (x, k±) 



m 



m 2 + k 2 \ 



(2.35) 



dx d 2 k i , , . , , , , . s M f ^ i xq± ■ k\ 



2(2tt) 3 

where k' ± = k± + xq± and 



-<p n >(x, k'^cp^x^kx) 



1 + 



m? + k 



2 / ' 



(2.36) 



M' 2 



m 2 + k' 2 m 2 + k'f 



+ 



(1-x) 

For charge form factor in the time-like region, 



x 



(2.37) 



d 2 k ± 
2(2tt) 



dx- 



Ml - r(m 2 + k\) 



4> n (x/r, k±) 



r(l — x) 



m 2 +k 2 . 



+ 



dx- 



x(l—x) 

;i_ r)(m2 + fc 2 ) _£(^) M 2 



V2Jm 2 + k] 



A, 



<j>tf((x - r)/(l - r), k±) 



x(l — r) 



9" 



(2.38) 



where r = p + /q + . One can find that, for the two terms in Eq. fl2.38j) , there remains one 
bound state vertex function which can not be substitued with the momentum distribution 
function shown in (|2.34| ). In fact, this is the nonvalence contribution arising from quark-pair 
creation flSyiSfl . We illustrate this situation in Fig.l (d), the process q — > qqq — ► qn is 



represented by the empty circle. Another (f>^(x, kj_) is not applicable for the empty circle 
because the light-front momentum of the quark is larger than that of the daughter pion 
(i.e. pi > p + ). This makes the task of calculating the nonvalence contribution considerably 
difficult. 
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For transition form factors in the space-like region, 

= / 9n V^Vt f dxd 2 k\ , . m 

F s (q 2 ) = 2^ / — — i ^(x, fc ± ) 



x 



2(2tt) 3 
1 



y^m 2 + /cj 



+ 



and 



(1 - x)^ + xM 2 - 2A; ± • g ± xgi + (1 - x)M 2 + 2A; ± • q ± 
W c rdxd 2 k 



m 



m 2 + k 2 



± 



x 



1 



+ 



1 



(1 - x){ql - q' 2 ) + xM 2 - 2k ± ■ q ± x{q\ - q' 2 ) + (1 - x)M 2 + 2A; ± • q ± 
For transition form factors in the time-like region, 

m 



3 J 2(2tt)3 ^-^2 



x ( J dx <f) n (x, k±) 



x 



+ 



+ 



k i 



(1 -y) (m 2 + kl) 
1 — x 



I d.r 0,(.r... , 
^ (1 -y)(m 2 + /c 2 ) 

/ 4>n{x, k±) 

Jl-v 



y[m 2 + k\ ) - (1 - x)[x - 1 + y)M 2 

ry X \ 
+ dx fa(x,k ± ) 2 7 r— 2 

jo y(m 2 + fcjj — x(y — x)M 2 J 



and 



F *f (9 ' 9 ) - 2_ ry 2727 



x 



3 7 2(2tt) 3 J m 2 + fc 
/•i-s/ 

y dx 07r(x, fcj 



+ [ dx tp n (x, A;_l) 

/ 0tt(x, fcjj 
ry 

/ dx 0tt(x, , 
Jo 



+ 
+ 



x(l -y) 


(1 


— y) 2 {m 2 + A; 2 ) — x(l — y — 


x)q' 2 




(l-x)(l-y) 




(1 


— y) 2 {m 2 + /c 2 ) — (1 — x)(x 


- y)q' 2 




(1 - x)y 





y 2 (m 2 + kj_) — (1 — x)(x — 1 + y)q 2 
xy \ 



y 2 (m 2 + kj_) — x(y — x)q 2 y ' 



where y = q + /p + and 



M 2 + q 2 - q' 2 - ^/(M 2 + q 2 -q' 2 ) 2 -4:M 2 q 2 
V = 2M 2 

7T 

We can easily check that F° (0) = FUO) = F wy {0). 
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III. NUMERICAL RESULTS AND DISCUSSIONS 



We now compare our results for form factors with the experimental data. Before doing 
this, we need to determine the parameters appearing in the wave function <f) n (x, k±). Of 
course, it is assumed that this wave function is universal and processes- independent. 

One wave function that has been often used in the literature for mesons is the Gaussian- 
type, 



^ ±)o=A ^ exp J_A_|, ,3,; 

where M = 4(7r/u; 2 ) 3//4 and k z is of the internal momentum k = (k±, k z ), defined through 

n ei-kz e 2 + k z 

1 — x — , x = , (3.2) 

ei + e 2 ex + e 2 



with ti = \/mf + k 2 . We then have 



xM m 2 , + k 2 



M = e 1+ e 2 , *,= -^— (3.3) 



and 



dk z = eie 2 . , 

dx x{1-x)Mq 1 ' } 

which is the Jacobian of transformation from (x, k±) to k. This wave function has been also 
used in many other studies of hadronic transitions. Besides this wave function which has the 
exponential form, we also consider the power-law type one 

(2 \ n 
WT^) ■ (3 - 5) 

where n is another parameter. 

For the pion case, we assume that the constituent mass of the u and d quarks is the 
same, i.e., m u = = m. Thus, there are two parameters m and u in the Gaussian-type 
wave function and an additional parameter n in the power-law type one. We will use three 
conditions to fix these parameters. The first one is the experimental value of the decay 
constant f n = 92.4MeV g| and (|2~35| ). The second one is the electromagnetic radius of the 



charged pion 

MM 2 ) 



(r ) 7r ~ -6- 



(3.6) 



dq 2 

where (r 2 )^. xp = 0.439 ± 0.03 fm 2 P|. The third one is the decay width T (tt° 77) K 



It is well known that F 7T7 (Q) can be determined from this decay width via 
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Thus leads to F£*{0) = 0.27±0.01 GeV _1 . Following from Eq. ( ggg ) we have the theoretical 
expression for F„- 7 (0) 



i^ 7 (0) 



N r f dx d 2 k 



2(2 



7T 



4> w (x, k±) 



m 2 + k 2 W 1 



x)M 2 



(3i 



Using these constraints, we can uniquely determine all the parameters in 4> w (x, k±) 
Here we list the fitted parameters of these two wave functions 4>c an d 4>n in Table 1. We 
show the x-dependent and /c 2 -dependent behaviors of these two types of wave functions in 
Fig. 2 and Fig. 3, respectively. Note that the location of maxima is quite different for <pG and 
4>n as shown in Fig. 2 and that the transverse momentum suppression of the exponential 
forms 4>g is more stronger than that of <pN (see Fig. 3). 



wave function 


m g (GeV) 


c^(GeV) 


A(MeV) 


(r 2 ) n (fm 2 ) 


i^ 7 (0) (GeV" 1 ) 




0.243 


0.328 


92.4 


0.434 


0.231 


(j) N {n = 1.7) 


0.192 


0.957 


92.4 


0.434 


0.272 



Table I. Parameters m q and w„ in wave functions 4>g and 



From Table I, we find that the value F 7r7 (0) = 0.27 GeV -1 cannot be reached by adjusting 
the parameters m q and u n in the wave function (j)Q. It is well known that 7r° — > 77 is not 
only the dominant decay mode (98.8%), but also relates to all the observed decay modes of 
the neutral pion (see below). Therefore the Gaussian wave function, in spite of being used 
in the literature for mesons, is not suitable to describe the decay processes of the neutral 
pion at least. Thus we will use the power-law type wave function <p(x, k±) N to calculate the 
form factors F n (Q 2 ) and F ny (Q 2 ). The form factors F n (q 2 ) and F^(q 2 ), can be calculated 
by Eq.( |2.36| ) and Eq. ( |2.39|) , and the results are plotted in Fig. 4, and Fig.5, respectively. 
They are in agreement with the experimental data. It appears that LFQM is valid up to the 
scale of order q 2 ~ —8 GeV 2 , while the constituent quark model is only applicable to the 
low-energy region. In the other approaches, |13| used the Gaussian type wave faction and 



the axial anomaly plus the PCAC relations to calculate the space-like form factors F^i^q 2 ) 
{q 2 ), their results are in agreement with data. The pQCD approach has also been 



and F^ 

used to calculated the form factors F n (q 

T2 



?s 

7T7 



(q 2 



in q < — 4 GeV region and F~ 
in —8 GeV 2 < q 2 < —0.5 GeV 2 region. In those regions, their results fit the data well. 

It is worth while to mention that tt° — > 77 amplitude comes entirely from the anomaly 
in the soft-pion limit ||32|| . The anomaly approach begins with this amplitude, and in 
terms of modified partially conserved axial- vector current (PCAC) hypothesis, it relates to a 
relevant axial- vector amplitude. In the soft-pion limit, the 7r° — > 77 amplitude is proportional 
to an anomaly term. The major difference between this approach and ours is that the former 
one doesn't assume an internal structure for pion, that is, the pion couples to the quarks 
simply via the vertex 75. The pion decay constant in the anomaly approach is regarded as 
an input value, in contrast with this paper, it could be theoretically calculated from the 
integral of a quark momentum distribution function ( [2.35| ). 
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Now we use the wave function (j) n (x, k±) ( |3.5| ) to calculate the transition form factors 

in time-like region < q 2 < 0.02 GeV 2 . The equation ( j2.41| ) can be used to get F^(q 2 ) 

with the results shown in Fig.6. This form factor is related to the differential decay rate of 
o 



k" — > 7 e + e by 

dTilT 



where 



ri 



7T 







7 e e 



77) dq 2 



A 3 ' 2 1,-^,0 G e tf) 



M\ 



(3.9) 



A(a, 6, c) = a 2 + b 2 + c 2 - 2(ab + bc + ca) 



and 



G e (q 2 



4M 2 
q 2 



1/2 



1 + 



2Mf 



9 



(3.10) 



(3.11) 



We take the integration of q 2 and obtain the branching ratio 



B, 



T(7l° 



7e r e 



7T 



77) 



1.244 x 10 



-2 



(3.12) 



This value agrees well with the data (1.25 ±0.04 ±0.01) x 10~ 2 [[34j and the weighted average 
of the two experimental results (1.213 ± 0.033) x 10 -2 |2"5| . In the approach of pure QED 
(neglecting strong- interaction effects), B ie + e - has also been calculated and is found to be 
1.185 x 10 -2 The approach of VMD [RBJ predicts the mass of the vector meson to be 
My = 0.821 ± 0.015GeV and also sets F^(q 2 ) -> F^(0) approximately. 

Another decay mode related to _F* 7 (g 2 ) is 7r° — > 7 + positronium. This process occurs 
in the Dalitz decay 7r° — > e + e~7 below slightly the lowest possible invariant mass of the 
lepton pair since the Coulomb binding energy is negative. Owing to 4m 2 -C M 2 , the ratio 
F* 7 (4m 2 )/F,|; 7 (0) ~ 1. Therefore the QCD correction to the branching ratio is very small 
and the branching ratio £> 7+P ositronmm becomes 



r(7r° — > 7 ± positronium) 



77) 



32vra 1 



4m 2 



E 

n=l 



4m 2 



0.6a 4 



1.7 x 10~ 9 (3.13) 



where ^ n Lj is the nonrelativistic bound state. The experimental value B^+ positT 
0.29) x 10~ 9 ||. 



±.84± 



Next, the equation of ( |2.42j ) can be used to calculate F^ 7 * (q 2 , q' 2 ) . This form factor is 
related to the differential decay rate of tt° — > e + e + e~e~ by 







e^e^e e 



r(7r° — > 77) <ig 2 c?g' 2 



q 2 q' 2 





K r (q 2 ,q' 2 ) 







x 



A 3/2 ( 1 



,/2 



V 'M 2 o'M 2 



G c (g 2 )G e (tf 



'2\ 



(3.14) 



After the integrations over g and q we obtain the branching ratio 
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B e + e + 



rfTT 



e 1 e 1 e e 



T(ir° — > 77) 



3.00 x 10" 



(3.15) 



which agrees well with the data (3.18 ± 0.30) x 10 5 [38|. The prediction of B, 
pure QED [|35j is 3.46 x 10" 5 . 

The tranistion form factor F nl * (q 2 , q' 2 ) is also related to another decay rate r(7r c 



e" l "e" l "e _ e~ 



111 



(for the lowest order in QED) by P5| 
r(7T° -> e+e") 



where 



TZ{p 2 



r(7r° — > 77) 



a 



71 



\_4M e 



2\ 1/2 



Ml 



2i 



[M 2 q 2 -(p-qf 



F^iqMp-q) 2 ) 



e 1 e 



(3.16) 



(3.17) 



J ~^q 2 (p- qf [{q - p e )2 - M 2 } F ny (0) 

It is well known that the unitarity bound for B(tt° — > e + e _ ) comes from the on-shell 77 
intermediate state. This state generates the model-independent imaginary part of 1Z [|(J 



Im ft(M2) 



(2 J 111 


"1-/3/ 




1 + Ar 



where Pl = l- 4M 2 /M 2 . Eq. 



gives rise to the unitary limit 



>B. 



unit 

e+e - 



4.75 x 10 



-8 



(3.18) 



(3.19) 



As for the real part of 71, we find, in Eq. ( |3.17| ), the integration is too difficult to do because 
that q 2 and (p — q) 2 don't have to larger than zero simultaneously. Thus Re 1Z can be 
estimated by applying Eq. ( |2.4U| ) in the soft pion limit, p —>■ 0. Then one obtains 

Re TZ(p 2 = 0) ~ -23.33. (3.20) 



In general, an once-subtracted disperson relation can be written for Re 1Z |4"T 

- 2 Im 7l{p' 2 ) 



Re 7l{p 2 ) = Re K{0) + — f°° dp' 2 - 

71 JO ( 



(p' 2 — p 2 )p 



12 ■ 



(3.21) 



We find that Re 1Z(M 2 ) = 8.93 and B e+e - = 5.98 x 10~ 8 . It is in agreement with the 
experimental data (6.3 ± 0.5) x 10~ 8 . The method of VMD yields the mass of the 
vector meson My ~ 0.77GeV and also obtains a value 6.41 x 10~ 8 which is consistent with 
experiment. 



IV. CONCLUSION 

In this work, the charge and transition form factors of the pion have been considered 
in LFQM. The parameters appearing in the light-front wave function are deterimined by 
experimental data of the decay constant, electromagnetic radius of the charged pion, and 
two-photon decay width of the neutral pion. Then we used these formulae and wave functions 
to calculate the charge and transition form factors of the pion and the branching ratios of all 
the observed decay modes of the neutral pion. The obtained results are all consistent with 
data. Table II listed the comparsion with experiments and other approaches. 
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Br 


exp't value [|25|] 


this work 


VMD 


pure QED 




(1.213 ±0.033) x 10~ 2 


1.244 x 10~ 2 


i.i85 x io~ 2 [g 


1.185 x 10~ 2 m 


^7+positronium 


(1.84 ±0.29) x 10" 9 


1.7 x 10~ 9 




1.7 x 10~ 9 


B e + e + e ~e~ 


(3.18 ±0.30) x 10" 5 


3.00 x 10~ 5 




3.46 x 10- 5 J35] 


Be+e- 


(6.3 ±0.5) x 10~ 8 


5.98 x 10~ 8 


6.41 x 10~ 8 [| 


> 4.75 x 10~ 8 



Table II. Branching ratios for the neutral pion in this work and in some other models. 



The time-like region of the charge form factor was considered to be applied to the decay 
process 7* — > tctt, the calculation was not finished because that we do not have a reliable 
estimate of the pair-creation effect. Contrary to the case of charge form factor, the time-like 
region of the transition form factors were studied completely. The major reason for this 
difference is that the pair-creation effect will exist when there are hadrons in the final state. 
Concerning the model-dependent part, although the choice of the light-front wave function 
is arbitrary, it is assumed to be processes-independent. Therefore it is desirable if more data 
in the low-energy region, where the constituent quark model is more applicable, are fitted. 
This is the reason why we use the power-law type wave function. The major differences 
between the Gaussian type and power-law one are, first, the locations of maxima is quite 
different for x; second, the transverse momentum suppression of the exponential function is 
stronger than that of the power-law (n — 1.7) function. We shall investigate the adaptability 
of this wave function and the treatment of Z graph for other hadrons and other processes in 
the near future. 
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FIGURE CAPTIONS 



Fig. 1 Diagrams for (a) a charged pion decay, (b) the scattering of one virtual photon 
and one pion, (c) one time-like photon and two pion coupling, (d) the illustration of the 
nonvalence contribution in (c), (e) the 7*7r° — > 7 vertex, and (f) the 7*71"° — > 7* vertex. 

Fig. 2 The x-dependent behavior of 4>g and <$n {n = 1.7) at kj_ = 0. 

Fig. 3 The /^-dependent behavior of 0g and 4>n (n = 1.7) at x = 0.19, where 4>g{x, 0) ~ 
</>jv(x, 0) at this location. 

Fig. 4 The charge form factor of the pion in small and large momentum transfers. Data are 
taken from p(J for small momentum transfers and [^] (empty circles), (filled triangles), 
and |23| (filled square) for large momentum transfers. 

Fig. 5 The one off-shell photon transition form factor of the pion. The dotted line is the 
limiting behavior 2f n (0.185 GeV) predicted by pQCD. Data are taken from ||24j| . 

Fig. 6 The y-dependent behavior of the function f(y), where f(y) = |F 7 J; 7 (g 2 )/F 7r7 (0)| 2 and 
y = q 2 /M%. The dotted line is the assumption of pure QED: f(y) = 1. 
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